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ABSTEACT 


A  llpeazised  theory  for  the  unstea^  notion  of  a 
partially  cavitated  flat  hydrofoil  in  two  dinenaions  is  carried 
out.  A  second  linearisation  procedure  is  dised,  based  on  ideas 
of  Tinman  and  Guerst,  to  obtain  the  unstea^  pressure  distribu¬ 
tion  around  the  tqrdrofoil  and  the  resulting  force  and  nonent, 
as  functions  of  the  cavitation  nuiri>er  and  Strouhal  number 
for  given  pitch  and/or  heave. 
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NOMENCLATURE 

(i  normalized  pressure  (or  acceleration  potential)  ■■ 

Poo"P 

-  j .  ,  vdiere  p  Is  the  actual  pressure  and  d  Is  the  density 

1  2 

a  normalized  cavity  pressure  (cavitation  number  *  U  ) 

G  Green's  function  for  regular  term  for  pressure 
I  oneohalf  the  length  of  the  hydrofoil 

I 

c  position  of  the  rear  end  of  the  cavity  (assuming  hydrofoil 

between  -i  and  i) 


i  velocity  potential 

U  free  stream  velocity 

V  oscillation  frequency 

0  ^  (reduced  frequency) 

M(t)  slope  of  hydrofoil 
B(t)  y  Intercept  of  hydrofoil 
f  force  perpendicular  to  hydrofoil 

m  moment  on  hydrofoil 


i 

1 

I 

\ 
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hydroelAstic  effects  of  unsteady  motion 

_ DF  A  HYDROFOIL _ 


lo 


1.  INTRODUCTION 

In  order  to  analyse  the  hydroelastic  effects  of 
unsteatty  motion  of  a  hydrofoil  it  is  necessary  to  determine 
the  pressure  distribution  on  the  hydrofoil.  The  problem  of 
particular  concern  here  is  that  of  a  partially  cavitated  hydro¬ 
foil,  experiencing  simple  harmonic  oscillations  in  vertical 
position  (heaving)  or  angle  of  attack  (pitching). 

In  order  to, solve  this  problem,  it  is  first 
necessary  to  determine  the  flow  pattern  around  the  oscillating 
hydrofoil.  Furthermore,  in  order  to  produce  a  solution  relatively 
quickly,  it  is  desirable  to  make  certain  simplifying  assumptions. 
First  a  two-dimensional  flow  pattern  will  be  assumed,  i.e.,  the 
hydrofoil  will  be  assumed  to  have  infinite  span.  Secondly,  the 
fluid  will  be  assumed  Infinite  in  all  directions.  Third,  thickness 
effects  and  Initially,  camber  effects  will  be  neglected.  Horeover, 
the*  angle  of  attack  will  be  assumed  sufficiently  small  that 
linearized  theory  may  be  used.  Finally,  the  wstea^  motion 
will  be  assumed  to  consist  of  small  oscillations  around  a  steady 
motion  so  that  a  "second  linearization"  may  be  made.  Therefore, 
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the  flow  pattern  Is  assumed  to  fluctuate  at  the  same  frequency  as 
the  hydrofoil  oscillation,  l.e. ,  higher  harmonics  are  neglected. 

Timman  [3]i  has  made  an  attempt  at  solving  this  problem. 

In  his  work  the  assumptions  made  to  complete  the  definition  of  the 
problem  have  consequences  idilch  are  physically  unacceptable. 

Guerst  [1]  has  formulated  the  problem  In  such  a  way  that  this  contra* 
diction  is  resolved.  However,  he  has  not  carried  through  the  problem 
to  any  extent.  In  our  treatment,  Timman 's  analytical  machine^  Is 
used  throughout,  but  the  physical  assumptions  used  to  determine 
the  unique  solution  are  those  of  Guerst.  Although  Guerst 's  set 
of  conditions  Is  not  beyond  qiiestlon  we  felt  It  is  the  best 
possible  in  the  ll^t  of  present  knowledge. 

To  solve  this  problem.  It  was  first  necessary  to  obtain 
a  steaify  state  solution  by  Timman 's  method.  This  was  done,  and 
the  result  compared  to  that  of  Guerst  [2].  Preliminary  calculations 
Indicate  complete  agreement,  except  for  discrepancies  vdilch  can 
be  attributed  to  numerical  Inaccuracies  (of  less  than  17»)  In 
Guerst 's  calculations.  Tinman's  approach,  as  developed  here, 
gives  much  simpler  closed-form  expressions  for  the  pressure  and 
velocity  potentials  than  Guerst  was  able  to  obtain. 

The  next  step  In  solving  this  problem  was  to  develop 
an  explicit  mathematical  representation  for  the  physical  assimptlons 
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made.  This  was  done,  and  resulted  In  six  simultaneous  linear 
equations,  Involving  six  unknowns.  The  coefficients  are  functions 
of  reduced  frequency  and  cavity  length.  Some  of  them  are  defined 
by  Integrals  which  must  be  evaluated  numerically.  To  do  this  a 
program  Is  being  written  for  the  IBM  7090,  which  will  compute 
these  Integrals  as  functions  of  the  two  parameters.  The  expressions 
for  pressure  and  velocity  are  simple  combinations  of  these 
coefficients  and  elementary  functions. (see  Section  4) 

Finally  the  lift  and  moment  of  the  hydrofoil  were 
calculated  as  linear  functions  of  the  above  mentioned  coefficients. 
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2,  GENERAL  STATEMENT  OF  PROBLEM 

Tinman  has  derived  the  general  expression  for  the  acceleration 
potential  ♦  to  be: 


♦  =  o  -  f(x*  ,t)G(x' ,o,x,y)dx'  +  Cj^#j^+C2* 

wetted 
surface 


(see  Figure  !)<. 

In  this  expression  c^  and  c^  are  functions  of  time  which  are 
determined^  together  with  the  position  of  the  rear  end  of  the  cavity, 
by  additional  conditions  given  by  the  physical  model  used. 

If  h(xs,t)  profile  of  hydrofoil 
then  w(x,t)  -  h  +Uh 

*  t 

and  f(x,t)  -  Wj.  +  Uw^ 

For  a  flat  plate ^  the  above  integration  may  be  done  explicitly.  (See 
Appendix  I) . 

Let  h  M(t)x+B(t) 

Then  w  =  M'x  fB'  +  MU 


and  f  -  M'*x  +  B’  +  2M'U 


Therefore 

♦  -  0  4  f  {B”  4  2M‘U)Fj^  +  +  C2^2 

M  (B’'  4  2M‘U)  ,  C2  are  functions  of  t  only. 


F^,  are  harmonic 
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P!{YSICAL  PLANF, 


5. 


6. 


functions  of  x  and  y  in  the  plane  slit  along  the  hydrofoil. 

Let  c  -  y  ~ 

where  z  •  x  -f  ly 

and  C  li) 

2i  ■  length  of  hydrofoil 
c  >i  x-coordlnate  of  rear  end  of  cavity 

_  „  /  i  -  c 

V  i  +  C 

Fk>  Figure  1)  are  harmonic  and  rational  functions  of  ^  and  t) 

only,  and  are  time  dependent  only  through  a. 

The  solution  of  the  general  problem  requires  specifying 
three  conditions  to  determine  the  three  parameters  c^,  C2>  and  a. 

There  are  four  conditions  which  may  be  used: 

Boundary  conditions  (general) 

(1)  ♦(«»)  ■  0 

(2)  >  w(x,t)  on  wetted  surface,  where  i  •  the  velocity  potential. 

(3)  h(c,t)  -  h(-f,t)  -  i  -  ^^)dx'  (closure) 

(4)  I  (x,0'^,t)  «  ly(x,0”,t),  X  >  i  (continuity  of  vertical 

^  ^  component  of  velocity) 
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Conditions  (1)  and  (2)  seem  quite  satisfactory  on  physical 
and  mathematical  grounds,  so  that  the  choice  lies  between  (3)  and  (4). 
For  the  steady  state  problem  (4)  is  automatically  satisfied  so  that 
(3)  is  the  only  condition  left.  However  in  the  unsteady  problem, 
the  physics  of  the  problem  seem  to  us  more  reasonable  if  (4)  is  used 
rather  than  (3).  This  is  the  view  of  Guerst,  and  others,  and  this 
approach  is  used  here. 

Using  condition  (1),  and  following  Tinman,  we  write 


-  •«  ks»l 


where 


c^  iV,  C3-/(B"  +  2M*U), 


1  X  4  X  4 

^  'k  v  ^  °k  ♦kx' 


o<»  k»l 


-  «•  k»l 


But 


^k  ^^k  1  ^^k  d»  x-x* 

^kx*  “  ^  ^  ■  U  "ST  d7  *  (“f  •  U“) 


•to  J 
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Therefore 


ks:!  “ 

The  harmonic  functions  are  listed  in  Appendix  VI,  Table  I 
where  those  with  subscripts  1  and  2  were  derived  by  Timman  and  those 
with  subscripts  3  and  4  are  calculated  in  Appendix  I. 
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3.  STEADY  STATE  SOLUTION 


Since  •  0,  we  have 


♦  -  0  +  +  C2#2 


*y  -  n  +  '2*'2> 


The  boundary  conditions  become  (1,  2,  3) 


(1) 


CiA 


-  C2A  +  *0  -  0 


c,B 


2  1" 

(2)  MU  +  """  +  CjB  >  0  ((  >  0  on  watted  surface) 

(3)  «Ve+l)  -  c,  /  (-**-»  -  ^  )»ti  +  c,  /  (t-B)dz 


where  A  and  B  are  defined  in  Appendix  VI. 
If  we  substitute  (2)  in  (3)  we  find 


(3‘) 


/ 


A. 


2^ 

-i  n  +r 


dx  c 


i. 


2  V  « 


-  0 


Since  X  ■  i  - 

u(r-Ki)^ 


and  T)  -  0  on  the  cavity  > 
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(3*)  is  replaced  by 


1  dx  /•**  »  j*  A 

'1  ?  3«  +  '2  {,  ‘  a?  •*«  -  ® 


or 


<=1  / 


iLls - H.  ^  „  f“  - t^idiSl.  .  d{ 


I  <1«  +  <=2  I 


[l+(t^-to)*) 


The  Integrals  are  evaluated  in  Appendix  III  and  we  obtain: 


(3*) 


iJJ+l  + 


2a) 


Ja^+1 


+  ( 


+1)  c,  -  0 


The  solution  of  (1),  (2),  and  (3*)  Is  as  follows: 

_  ^  ^  2U^ 

Let  p  - — 


Consider  the  eqiiatlon  y  -  y  +  p  »  0 
Then 


Let 


o  -  i  4  -  y) 


-P+l  /^  - 


2  .  4 
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n. 


Then 


T  -  + 


3X 


I 

173 


The  only  allowed  roote  axe  those  which  give  positive  y.  In  general 
there  are  two  solutions  y^,  ^2 

^  ^  ^  ^  solution  which  is  physically  neanlngful 

\/3 

(see  [21)  is  given  by  y^. 


Special  Cases: 


1) 


2  4 

p  •  Yf  (eppetcncly  the  nsxlnun  angle  of  attack  for  which 
linearised  theory  will  hold.) 
y  «  (2  positive  roots  identical) 


c  «  — |— 

«  .75  (cf.  Geurst  value  of  .7485) 

2)  p  «  0  (no  angle  of  attack) 
i  fiL  £ 

0  *i  (no  cavity,  physically  real  solution) 

1  0  i  (cavity  full  length;  anomalous) 
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Since  y  is  only  an  Intexmediate  variable,  it  is  of  interest 
to  express  everything  in  terms  of  o. 

By  solving  ^  -  y)  t  we  obtain 


Since  p  >  y 


«  -oB(l-fa  -  o) 

«  oB(a/oi^+l  +  1  -  o^) 
p  ■  2aiJa^+l  -  o)^ 

Therefore 

Cf  -  (o^+l)  (o  +  Vo^+l) 

e,  .  (3,^  3 

A  o 
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4o  UNSTEADY  CASE 

The  general  expression  for  the  velocity  Is  given  by 

U»y  =  2  |cj^(t)^j^(x,y,o(t)) 

Cj^(t-2^)4(x',y,ci(t-3^))dx' ] 

-ip 

Assume  the  time  dependence  Is  given  by  a  small  simple  harmonic  oscillation 
around  the  steady  state 


a  ®  oIq  ®1  cos(vt-0^) 

‘k  -  +  ^ki  ““(''C-Sit) 

where  a^.  Cj^q,  C2q  are  the  steady  state  values. 

Let  M(t)  =  +  Mj  cos  vt, 

B(t)  «  cos(i/t-0g) 

where  M  ,  M, ,  j  6.  are  all  given  quantities, 
o  X  JL  n 

Then  M»  (t)  =  -  sin  vt 
ri"(t)  =  -  v^Mj^  cos  vt 
B”(t)  =  -  cos(vt-0g) 

=.  -  vi(vBj^  cos  (vt-0g)  +  2Mj^U  sin  vt) 

=  -  v>?(vBj^  cos0g  cosvt  +  (2Mj^U +vBj^  sin  Sg)  sin  vt) 

- - - - - -  ?M,U 

=  -  vi  )f(vBj^)2  +  4(Mj^U)^  +  AvBj^Mj^U  sin  0g  co8(vt-tan  ^(tan 
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Cjj^  .  -  Vi  /(vBj)^  +  4(M^0)^  +  4vB^Mj^0  tin  »g 

.1  “l® 

9-  •  tan  *"  (tan  9-  +  '  '  ■  ) 

^  *  vBj^coa  6^ 

COi  vt 


Tha  problM  la  to  find  Cy^^,  c^y^,  9j^»  S* 

Wa  now  eavxy  out  tha  aaeond  linaavlsatlont 

(Aovav  a  ayabol  danotaa  tha  paxtuzbatlon  In  that  quantity  (Ivan  by  tha 
aynbol  without  tha  A) 

Ck(t)t|j(».y,a(0)  -  Cko^k^**y*®p>  CkiCoa(vt-9k)tk(**y»®o^ 

o  . 
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^  *y  ^  (x, y, 0^)03^  cos ^t  -  0^) 

■y  k«l  ^  .  o 

■*■  u  I'^ki  V^k^*'  »y»®o^ 


-iP  I 


“kO  U'.y,ao)c.jSln(v(t-S^)-«^)j  dx'l 


The  above  Integrations  are  carried  out  explicitly  for  the  cases 
where  y  ■  O’*"  or  y  =  0",  thus  giving  the  vertical  component  of  the  pertur¬ 
bation  velocity  along  the  hydrofoil.  The  steady  state  term  has  been 
derived  previously.  Similarly  the  folkxfling  equation  gives  the  pertur¬ 
bation  term  for  the  pressure  distribution  along  the  hydrofoil  when  y  k  O'** 
or  y  «  0“ : 

A  ^ 

k-l  o 

of  ^2  given  in  Appendix.  71,  Table  I, 

o 

of  are  not  needed  since  c^q-c^q-O. 

o 

The  boundary  conditions  (1) ,  (2) ,  (A)  of  Section  2  are  explicitly  given 
as  follows : 
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Condition  (1) 


4 


£  (C|^j^cos(vt-8^).J^(-)  +  3^  t,jcos(ve-e^))-0 

k.“l,  o 

where  the  values  at  «  of  and  t —  are  given  in  Appendix  VI,  Table  II 

K  O  Ot 

Condition  (2) 

w  M»x  f  B'  +  U(fM-M  ) 

'■  xvMj^sir-  s"  -  .  Bsin  (i/L-bg)  4  UM^c  vt 

-  -  xvM^  sinvt  +  (UM^ -}  vB^  sinfcg)  cosi/t  ^  vB^^cosS^sin  vt 

^41  U  ^41 

— 77  (x  sin  >t  4  “  cos  7-'t  )  4  “““•  sir(vt-@^) 

^,^.2  V  vi  3' 

Therefore  Condition  (2)  becomes 


Uc 


Uc. 


C|  sin  .-t  -f  cos  vt)  f  “““  sin(vt-03) 


4 

2 

k..3 


X 


^  ^'kl  ^  /  V'kCx  .0“,o^)sin(v(t- ^^^)-0jdx' 

-At 


2  cos  ivt-ti^)i>^{O.J)\aJ-i  C|,n  (OsO",a,^)aTCOsO^t-0^) 


k=.d. 


O"  kO  da 


o'“'r 


I  r  [c 

u  J  L‘- 

*00 

hij 


4-  C, 


kO  da 


—  (x*  ,0  gO  sin  «'v't  -  )  dx'] 

8  j  O'  1  •  U  a'' 
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■  ^  Cj^llJ'l^(O»O",aQ)[co8(vt-0j^)-co8(v(t— 

lc*l 


'itO  5?  (0,0',a^)o^[cos(^e^)-co»Kl;  -^)-93)] 


[♦ 


,  -r;^  axe  con8taet8  along  the  wetted  eurface  for  k«l,2j. 


k*  da 

Thl8  may  be  elxapllfied  to 


4  *11. 

2  Cj^l(^lj(x,O‘,Oo)co8(vt-0k) /  8tn(v(t  ,0’,o^)dx*) 

lcv3 


•«0 


1  , 

j  dx' 


k  .  ,  «  X  xtx' 


+  “I'TcO  55"  (x'.0,a^)8in(v(t  — 
o 


-  x+1  V  ''^k 

+  Ckl^k(0»0  »<»^,)co8(v(t  -n|r^“V‘*'®i*TcO  do”  ,0^)c68(v(t- 

O  ' 


Uc 


Uc, 


^  sin  vt  +  ^  cos  vt)  +  -r;^  sin  (vt-03) 
-i 


•^k,c  “  &  ^  ^k^*' 

-«9 

4,0-  4.,  <*®*^*<*  5^ 
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These  Integrals  must  be  evaluated  nximerlcally . 


The  f,  and  ^ —  needed  are  listed  in  Appendix  VI,  Table  V  (in  the  form 

§ci~*  non-constant  parts  of  V'jj  30”*  i*’  Appendix  I) 


Let 


Then 


p  -  v(t-  ^  ), 


w  =  /  sin(p  +  ^  ,o~  ,a^)dy.* 

"  i 


1  ^  VX* 


-  j  /  sin(p  +  -  e3)(x»+i)dxV 

i 


-2cos(p^.J^-9 cos(p  4.  Jf-  -  93)dx' 

-£  -i 


X 


/ 

-i 


8in<p  +  ^  -  03)V^3(x'  ,o,a^)dx* 


2 

-  -  ^  V'3U^o"»a  )co8(vt-02)+  ^(8in(vt-02)-8l-n(v(t-  ^^>-^3)) 

V  I 
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IS 


Lee  T  =  8ln(p  +  (x*  ,0-,o_j)dx' 

-  -a  2  cos<p  +  J^)[(i^)2  .  (2^)  j* 

"  i 

^  2  ;  e„3(p  .  1]  ax. 

-i  2i^ 

-  -  2  2  cos  (p  +  :^)  [(2^)2  .  (2^)J 

.22  /%os(p  .  22^)  ^  ax.  .  T,  .  X, 

^2  '  ^  ]*,  -  /  >1"<P  +  t)  *.  I 

^  -J  -J  U  J 

=  ^  [x  sln(p  .  /  si„(p  -  ^)  +  2  (eoaCp  .^)  .  eo.(p  - 

Note  that  =  -  2  CO.  (p  .  ^) 

Therefore 


X 

sin(p  +  ^)  (x»,o",a^)dx» 

“  ■  V  1*'4<*^»0",a^)fco8(p  +  ~)  -  cos(p  -^)) 
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20. 


=  -  ^  V'4(x,o",a^)  cos  vt  +  (j  sin  vt  +  :^  cos  vt) 

Jv 

2 

+  (sin(p  cos(p  -^))  +  “  %  (-i.o’,ao)cos(p  -  ^)) 

V  i 

Condition  (2)  then  becomes 
where  s  =  t  - 

“l'^kot^‘^k,s+^  (-i.o‘,a^))cos(vs-e^)+ J*^^sin(vs-0^)] 

■  ""ai  n  ®^”(''®>  •  (^)^cos(vs)) 

Condition  (4) 


4  i 


®  ^  ^  Cj^j^sin(v(s -0j^)  [v^j^(x»,o'^,a^)-V'j^(x',o",a^)] 

k=l  -i 

+  Cj^^aj^sin(v(s  +  i|jt^)-0^)[^  ^  (x' ,o”.a^)]  dx' 


Integration  along  the  wetted  surface  in  the  z  plane  is  equiva¬ 
lent  to  integrating  along  the  q  axis  in  the  ^  plane.  Similarly  Integrating 
along  the  cavity  sur!face  in  the  z  plane  is  equivalent  to  integrating 
along  the  ^  axis  in  the  ^  plane.  Therefore  (4)  can  be  stated  as  follows: 


TECHNICAL  RESEARCH  GROUP 


21. 


4  p 


x*+i. 


2  /  Cj^^8ln(v(8  +  ^^)-0j^)V'j^(o,n)+Cj^Qaj^8in(v(s  + 

1-00  o' 


x'+i. 


^4, 


+  /  C|^i8in(v(8+ ^^^)-flj^)^j^(5,o)+Cj^QOj^8in(v(8 +^^^)-9^)^(4,o)-j^4-0 
o  o  ’ 

Let  •  non-con8tant  j^art  of  Since  the  above  equation  la 
satlafled  by  the  conatant  part,  It  la  aufflclent  to  conalder  only  the 
non-conatant  parta,  which  are  given  In  Appendix  VI,  Tablea  III  and  IV,  for  the 
wetted  cavity  aurfacea  reapectlvely. 

On  Inapectlon  of  theae  tablea  we  aee  that  along  the  wetted 
aurface  and  all  equal  0  for  k>l  and  2.  Furthermore  for  k«3  and  4, 
on  the  wetted  aurface  are  both  polynomlala  In  x,  while  on  the  cavity 
aurface  are  equal  ^  quantity  defined  In  Appendix  I,  plua  Che 

aame  polynomials. 

Therefore  condition  (4)  can  be  almpllfled  to  the  following: 


°  ?  J'"  8ln(vs-9^t  d 

1-00  ^ 


Vjt 


where  0  a  (reduced. frequency) 


P  .  I  +  («^  +  o)^  , 


and  la  replaced  by  for  k-3  and  4. 


TECHNICAL  RESEARCH  CROUP 


22. 


Let  I 


k,c 


I 


J 

•i«Q 


CO  s  dj 

♦^sin(f)  lii^de 


where  Is  used  instead  of  for  k>3  and  4. 

*  dV'u 

Let  I,  be  the  corresponding  integrals  involving 

o 

These  integrals  are  reduced  to  infinite  series  in  Appendix  IV. 
Then  condition  (4)  can  be  stated 


*  c 


Since  the  conditions  nnist  hold  for  all  time,  each  condition 
is  expressed  as  two  linear  equations,  given  by  the  coefficients  of  the 
cosine  and  sine  of  vt  or  vs,  respectively. 


"'k.c  ■  ‘kl'®*  ®k 

'*k,8  -  ®kl**“  *k. 

a  a  a.,  cos  9 
c  la 

“8  '  “l®‘"  «o 
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23. 


Condition  (1)  becomes 


1 


f  |Vc  +  ‘^ko  -ss;  “c| 

f  |\,s  +  ‘=kO  -S^  <“>  “b  I 


Condition  (2)  becomes 


4 

2 

1 


\,c  \,c  \,s^‘^k,s  »“o^)  ■'■  ^kO  ‘^k,c  “c 

*  ^\r  -  T 

+  ''kO  ^‘^k,s  +  (-i.o‘,aQ))ag  / 

fS  ^ 


Vi  *^3,0  ■*■  Vi  ^vi^  ‘^4,s  ”  ° 


4 

2 

1 


“‘^kjS  ‘^k,c  ‘^,c^‘^k,s  ■'■  ’“o^^’^^kO  ‘^k,c  ®s 

V-  -  *) 

■*■  ''ko<‘^k,s  ■'■  •ss;  »%»“c  j 

+  ‘^3,s  “  H  "^4,8  +  *^4,0  “  ° 


Condition  (4)  becomes 


f  |^,c^k,c  ■*■  \s^k,s  ■*■  ®k0^k,c“c  ‘^k0^k,8®sj' 

f  I’Vs^c  ■’■  Vc^k,s  ■  ^koVc“s  ■'■  ^koVs^cj" 


-  0 
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24. 


After  solving  the  above  set  of  equations  fox  d.  ,  d.  ,  d_  , 
d  >  a  ,  and  o  ,  we  then  can  explicitly  write  the  unstea(fy  pressure 

2  g  8  C  8 

term  ae: 

*  ■  [Ji  VA  “c^'kO  ^)J<=o»v‘ 

From  this  the  unsteady  lift  and  moment  can  be  obtained. 

1  . 

These  are  explicitly  developed  In  Section  5. 

Since  the  closure  condition  was  dropped  from  the  calculation 
It  Is  of  Interest  to  obtain  the  cavity  shape.  To  obtain  the  horizontal 
component  of  velocity  along  the  hydrofoil,  needed  for  such  a  calculation, 
we  would  have  to  evaluate  certain  Integrals  which  are  not  otherwise 
needed.  Therefore  the  solution  of  this  problem  Is  postponed. 
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25. 


5.  THE  FORCE  AND  MOMENT  ON  THE  HYDROFOIL  (UNSTEADY  CASE) 


Let  g(x,t)  «  ♦(xjO.t)  -  ♦(x,0‘*‘,t) .  Since  ♦  -  (p^o-p)/*!, 
where  d  is  the  density,  the  pressure  difference  at  each  point  x  of 
the  hydrofoil  at  time  t  is  given  by  -dg(x,t).  Thus  Hhe  force  f  normal 
to  the  foil  and  the  moment  m  (around  the  leading  edge)  are  given  by 
(see  figure  2): 


f  m  -d  f  g(x,t)dx 
-i 
i 

m  -  -d  p  (x  +  i)g(x,t)dx. 

l^e  lift  «  f  cos^,  where  P  is  the  angle  of  attack  (given  by  cosp 
The  moment  m'  around  any  other  point  Xq  is  given  by: 


m'  (x^)  »  m  *  (x^  +|)f. 


.  A 

Since  ♦  =  Sq  +  ♦, 

A 

and  since  ^  can  be  expressed  by 

A 

4  =  t^cosvt  +  ♦^sinvt 
c  s 

where  4 

♦c  -  ^,c%  “c^'^kO 

4  d4, 

*a  *  \8*k  +  ®s^®kO 

it  follo\«s  that  we  can  similarly  resolve  g(x,t)  and  the  resultant  f  and  m: 
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pivot  about 
which  the 
Moaent  is  taken 


Force  (£)  and  aonent  (a)  are  positive  in  direction  of  arrow 


Figure  2 


27. 


g(x,t)  =  g„(x)  +  g  (x)cosvt  +  g  (x)8lnvt 

U  w  9 


f  ■  f  +  f  cosvt  +  f  slnvt 

O  C  8 

m  »  m  +  m  cosvt  +  m  slnvt 
o  c  s 


For  problems  of  heave  alone  or  pitch  alone,  £  Is  the  magnitude 

c 

of  the  unsteady  force  In  phase  with  the  displacement  (assume  for  heave 
alone,  d-  -  0),  while  -£  Is  the  magnitude  of  the  component  of  force  In 
phase  with  the  velocity.  Similarly  m^  and  -m^  are  the  components  of  the 
moment  In  phase  with  the  displacement  and  velocity.  We  have 


«o  -  “k0<*k  •  ’k^ 


*c  ■  \,c^*k  ■  *k^  “c'^kO  ^ 


_ k  _ kv 

o  o 


8s  -  \,s^*k  ‘  ‘^k^  ®s\0 

Let  Rk  -  2  /  (♦■  - 


£ 


Sk.2  /  -  ♦;j)dx 

*  * 

and  let  R^,  be  defined  similarly  using  Instead  of  The  Integrals 

o 

above  can  be  considered  as  llife  Integrals  (In  the  (x,y)  plane^}  ovet  the 


surface  of  the  hydrofoil  In  a  counter  clockwise  direction.  This  Is 
equivalent  to  Integrating  (In  the  (^,t))  plane  -  see  Fl^re  1)  from  -oeto  0 
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along  the  tj  axis  and  then  Integrating  from  0  to  c#  along  the  ^  axis.  Since 
all  and  ^  are  zero  along  the  §  axis,  Sj^  and  Sj^  can  be 
expressed  by 

Rk  -  2  /  \iOM 

*40  * 

•CO 

h  -  1  (*+ ^  an 
*  .  *  *" 

^  ***  ^**'*"*^  replacing  by  ^  in  the  above.  Note  that  all 

^k*  SaT*  while  x+i  is  even  as  functions  of  tj.  For 

0  ' 

purposes  of  integration: 


X  +  i  ■ 


XrT  -  +  1 

-8iii(Ti^-  a^) 


((TJ  -O^)^  +  IV 


The  integrations  are  carried  out  in  Appendix  V. 

I-etR  .  JI  c  V; 
k-1 


Z  c.„S 


kO“k 


TECHNICAL  RESEARCH  GROUP 


29. 


Then 


0  ^ 


k=l 


kO“k 


f  = 


=  T 


-d 


T<W  *  '*k.A> 


and  in  ,  m  ,  and  in  defined  similarly  using  S,  and  S  Instead  of 
o  c  s  K  a 

R,  and  R  .  f  ,  m  ,  R  ,  and  S  are  calculated  In  Appendix  V.  The  steady 
k  a  o  o  a  a 

state  terms  are: 


f^  =  -irdiU 


(1  +  \ll+  ) 


“  * - (Va^+l  4-  a)^) 

°  Aa(o^+l)  ^ 
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6.  NUMERICAL  ANALYSIS  FOR  UNSTEADY  SOLUTION 


The  "l"  integrals  may  be  evaluated  using  the  results  of 
Appendix  IV,  while  the  "j”  integrals  require  direct  numerical 
integration,  using  Euler’s  method"  for  increasing  the  rate  of  con¬ 
vergence  of  alternating  series.  Given  these  integrals  in  terms  of 
and  Q,  one  problem  is  reduced  to  solving  the  system: 


Qv  =  q 

where 


and  and  Q2  are  each  3x3  matrices: 
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where 


31. 


it 


-  4-  z  =^o  ^ 

a  ^  kO  da 


IfC”) 


^^(<») 


(oo) 


U  Mo  1 


^  - - Is- 


'kO 


1  ^  * 

U^M  1  hO  k,c 

o 


Similarly  and  1^^^  are  defined  in  terms  of  J*^^, 

^k,s*  ^  <‘^1,8*  4,8»  “s»  ‘^l,c»  4,c*  “c> 

vector  q  is  given  by: 


‘^k,c<‘^k,c  +  +  *^,8(^,8  +  ^  ‘ 

^,c\,c  '*'  8^k,s 

'C.c^k^") 

‘^,c(•Jk,s  +  ^  -  ^k("»-4,s<\.c  + 


4 

z 

k^3 


^,c^k,s  "  ‘HcjS^kjC 
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where 


Since  any  solution  to  a  problem  of  pitching  and  heaving  can  be 
expressed  as  a  conibination  (taking  into  account  phase)  of  pitch 
alone  and  heave  alone,  it  is  necessary  to  find  only  the  solutions 
for  unit  amplitude  pitch  alone  and  unit  amplitude  heave  alone. 

Let  qp  and  be  the  unit  pitch  and  heave  vectors  respectively, 
given  by  for  pitch  alone  and  for  heave  alone. 

For  pitch  alone: 


Cjj^  "  *  2v£M.jU 


03  -W2 

C41  -  - 
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Therefore  (for 


<13,0  -  0 


'3,s 


*14,0  - 

.*  - 

\s  ■  ® 


2a 

2 


s 


+  5)  +  “<J3.s 

1 

ofl 

/-V 

8 

^\c  •* 

1 

HJ- 

1 

8 

Km^ 

1 

2^3.0 

"^14, s  - 

2ain  ^ 

3,c 

For  heave  aloiie: 


'31 


-  u^) 


O3  -  0 


C41  -  0 
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Therefore  (for  B-j^=  /M^) 


,*  2 
3,c 

A*  _  A* 

S,s  ^4,c 


d,  =0 
4,s 


= 


Let 


.q\ 


3,c 


"^<^3,3  - 


3,s 


''s  =  K,S'  4.S’  “3> 
''c  =  (“I.c-  ■*2.C-  “c> 


and  let  v  ,  v  be  the  unit  solutions  for  pitch  alone  and 
s,p  c,p 

V  ,  ,  V  ,  the  unit  solutions  for  heave  alone.  Let  v  .  ,  v  , 
s,h*  c,h  s,h*  Q.,h 

be  the  unit  solutions  for  heave  alone  and  0^  =  Tr/2.  Then 


,  =  -V  ,  and  V  ,  =  v  ,  .  The  unit  solution  for  heave  alone 
c,h  s,h  5,h  c,h 


at  an  arbitrary  phase  0  is  given  by 

s 
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.cos  0_  +  .sin 
s  Syti  B  c,n 


V  V  .cos5„  -  V  .  sin0^ 

c  c,h  B  s,h  B 


Therefore  for  unsteady  pitch  «  M^cosvt  and  xinsteady  heave  * 

M  B 

''s  -  M-  ''s.p  +  rr 

o  o  * 

Ml  Bi 

V  «  —  V  +  TT-T  (v  .  COS0_  -  V  .  sind_) 
c  M  c,p  M  i  '  c,h  B  s,h  B' 


Bicos(vt-0g) 
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APPENDIX  I  -  fy  ^3,  04 


Calculation  of  ^3,  4*3,  ^'4*  ^4 
A.  General 


-4:5  /  x'''^0(C,C’)dx 

2Tri  wetted 
surface 


Where  G  ■  Re(in($  )) 

_ C'-C 

■/ \/M  ■  “ 

Along  wetted  surface  Is  imaginary 

?■  -  ^flW- 

and  X  =  - - r  -  i 

1  +  (u^-a)^ 

Therefore  G  =  in(^'*'H  •  ,  since  argument  of  in  is  a  positive  real 

u+iC 

nvimber . 

As  a  result,  G  =  F(5)  +  F(5) 
where  F(0  = 

Since  as  a  function  of  the  are  the  imaginary  parts  of  an  analytic 
function  of  where  are  the  real  parts,  each  can  be  expressed  as: 
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\  "  2i 


Let  H(0  =  2i?(0 

Then  G(0  »  ^  (H(0  -  H(C)). 

Therefore  f'  *  — f  H(5,lu)  dx 

2wi  wetted 


surface 


-i 


wi 


k-2 


where  ~  being  determined  only  up  to  an  additive  constant, 

so  that  -  -«’)  -  0. 

The  wetted  surface  of  the  hydrofoil  corresponds  to  the  interval 

’'OjOo)  of  u. 

Therefore  /“  ^  /n(^)  du 

trl  o  ^ 


,  -1  fn(H4 

(k-2)Tri'^  ^ 


(k-2)Fi 


i  *  k-2  ,  .vk-2,  1 

^  I  ^  -  (*^)  ( 


u+i?  u-iC 


)  du 


At  u  -  00,  X  «  -/,  »  0 


at  u  .  0,  X  ---  i 


l-g 

Ita' 


-  ’Ti 


Since  is  only  known  up  to  an  additive  constant,  the  first  term  for  f|^ 

Fay  be  dropped ^ 
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Let  be  the  non*constant  part  of 


Since 


,  the  Integrand  Is  even  and 


fj^  (=  f,'  -  above  constant)  can  be  represented  by: 


I  (- 


k-2 


■k-2(k-2),r 


7  -  1)^  -(-1) 


-k-2 


^u+iC  ”  u-i^^ 


du 


In  particular! 


'•u+iC  u-iC 


)  du 


f/.  = 


i  r 

z  I  !< 


TT 


oo  L  l+(U^-a)^ 


l+(u^-a)^J 


'u-iC  u-iC 


)  du 


Let  t  =  A+Bi,  A  and  B  positive. 


where  t  =  a  +  i 

7  2  2  7  2  -2 

Then  1  f  (u  -o)  *  (u  -t^)(u  -t  ) 

Since  ^  *•  where  ^  >  0. 

Therefore  for  purposes  of  determining  residues  to  evaluate  f^  and  f^  by 
contour  integration,  the  upper  half  plane  poles  are  t,  -t,  and  i^. 


Since 


'  (■  ' 


t  1  2.,  2  -2.  “  ,.2  -2  '  2  .  2  2  -2 

(u  -t  )(u  -t  )  t  -t  u  -t  u  -t 


) 


=  ~(  ^ 
04  V 


2i  ^  2  2  2  -2^ 

u  -t  u  -t 
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1  oo  1 

^  oZ  f  (“5 


3  ■  2'  '»V-t2  u2.e2'  'U4it  U-IC 


)  du 


f/.  =  -  -  / 

L  TT  J 


•oQ  L 


i  (-i- 

i*  ^  2  2 


u  -t  u  -t 


1  >2 

2  r2-'  .  2  ^2w  2  r2. 


(u  -t")(u  -OJ 


Wi?  u-i^ 


)  du 


-i 


1  x2 


‘  ^  ^-^u2-l:2  u2-c2’  ^u+lt  ■  u-lC>  ^ 


In  addition  to  the  general  expressions,  the  following  special 
cases  are  desired 

1)  at  X  =  <30  (or  ^=B  and  ti=-A)  ,  ♦  and 

since  the  desired  ^  is  0  at  this  point  and  the  function  must  be  adjusted 
accordingly 

2)  at  ^  =  0,  (the  wetted  surface)  ♦  and  ^ 

3)  at  r)  »  0,  (the  cavity  surface),  f,  and  separated  Into  odd  and 

even  functions  of 
2  2 

4)  ■'\^  =  ^  40 (for  J,  and  J,  Integration)  ^ 

cC  •  C  K  A  8 


B 


Calculation 


oo 

/ 

—  PC 

Zt'u-t 

u+t' 

-  ~(— 
2t  u-t 

1  ! 

^  2  ] 

1  1 

f  -2  ^ 

2t  1 

[c+icj 

t  2t 

1 ic-t  J  . 

u+t  ' 


U-1^ 


•t)  du 


t(t-lO 


tCttijll 


tt  I  (t+iO(t-ic)  (t-iO(t+ic) 


TECHNICAL  RESEARCH  GROUP 


40. 


t+iC  I  ^  -  |A+Bl--n+U|  ^  =  (A-n)^  +  (B+0^ 

t+ic|^  =  |A+Bl+n+t||^  -  (A+t))^  +  (B+0^ 

t(t+lC)  -  (A+B1)(A+ti+1(B+0)  -  A(A+Ti)-B(B+p+i(2AB+A4+Bn) 

■  a  +  ATi-B^+i(l+A|+Bii) 

t(t-10  -  (A-Bt)(A-n-i(B+0)  -  o  -  An-Be-i(l+Ae-Bn) 

1 1  *  ^^+1 


Therefore  f. 


Therefore 


^3 


±  At-  ^  ..  .1.  ii  A^^Bn  ^ 

+1  (A-ti)"^  +  (B+0^  (A+n;  +  (B+0^ 


y/o^ 


o-An-Bg 


_ o  .+  An  -  BB 

(A-n)  ^  +  (B+0  ^  (A+n)^  +  (B+e)  ^ 


JLjlM  ( _  1 _ ^ _ 1 

^  (A-n)  +  (B+e)^  (A+n)^  +  (B+0‘ 


■) 


btl 


g-Bg 


V(A-n)^  +  (B+0^  (A+n)  +  (B+0' 

) 


•) 


\(A-n)''  +  (B+O'^  (A+n)^  +  (b+0' 


( - 2^ 

V (A-n)^  + 


(B+0^  (A+n)^  +  (B+e)^ 


) 
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Special  cases 


1)  V'' 


1 _ +  + 

2.„2.  ^ 


2/Ai  4B 

1  1 


2  \/a^+l 


.4(A^+B^)  4B^ 


2(ci^+l) 


2  \/a^+l(  \/a^+l  -  a) 


■+1  -  a 


♦  .  = 


2(a^+l)(  ^a^+l  -  a) 
a^+2  +  a  \/a^+l 

2(a^+l) 

■3a-  L^+1 

( 

1 

2  \Jahl  ' 

^4 

a  2 

4^ 

/a^+1 

a  +  ' 

L 

2(a^l)  '  2  ^ 

TT 

1  +1 

a  +  7o^+l 

+ 


+1  +a 


2; 


4B' 


2p»I 


2)  >3 


2  (0*^+1) 
1  f 


4\KI 


4b‘^ 


_ 

2  I  ITT  1  :x 


\b^  +  (A-n)^  +  (A+n)^;  0^+1  VB^+(A-n) 


2(B^+A^+n^)  -  4ABn^ 

^yKI((A^+B^+Ti^)^  - 
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42. 
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where  f*  .  4  2,  - 

A)Let  D=  ((A-Tl)^+(B+0^)((A+n)^+(B+0^)  -(A^-n^)^+(B+e)^+2(B+|)^(AW) 

2  2 
Substitute  ^  -fa  fox  r\ 

D  -  (B+0^  +  2(B+4)^(A^+4^-»a)  +  (A^  -  o  - 

«  +4B^5  +6B^4^  +  4Be^+|^  +  2(B^  +2B5  +  -f  a  +  A^) 

4  2  4  2  2  2  2 

+  A^  +  a  +  I  +  2ar  -  2A  o  -  2k  ^ 

-  +  1 4  -  2o^  +  4Be(a  +  Jahl)  4  |^(6B^  4  2o4  2A^4  2B^4  2o  -  2A^) 

2  4 

4  8Br  +  45 

-  1  4  4B§(a  4p^)(4  §^(4  ^0^41)  4  '8B|^  4  4^^ 

~  2  ((14A5)(>/^  4  25^4  0  4  2Be)  -  (C^^g)) 

^3  _  _ 

^l-K*^(l4-4B5(a  4/0^4!)  4  45^(^^)4885^445^) 

,  2(25^4N/o^4  24(^0^41 -0)5  4A5(/o^4l  40  4  25^)) 

(l4  4B5(ci  4  ^0^41)  4  45^/0^41  4  8B5^  +  45^ 

=  2(26^4\/£?U  4  A5(3  yo^4i  -  o)  4  2A5^) 

yi^Cl  4  4B5(o  4/0^41)  4  45^  /o^4l  4  885^4  45^) 
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44. 


^'3 


flKt 


2  i 


(J^+l  +  2^  +  A£(3Va^ -  o  +  2ih)  Ja^+l  +  ^Hict-¥jahl+2i^) 
(1  +  vKl  +  4^^)^  -  16B^e^(o  +  \/?+l  +  2i^)^ 


=  1  -  4a^^  -  44^ 


7L  +  t 


(1  +  44^  /a^+1  +  44^)^  -  (44^  +  4a4^-l+4(  -  a)4^  +  2)^ 

-  (44^  +  4o4^-l)^  +  4(44^+ 4a4^  -  1)(2(^^^  -  o)4^  +  D 
+  4(2(  ^0^+1  -  a)4^  +  1)^ 

=  (44^  +  4a4^  -  1)^  +  4(2(vKI  -  a)  4^+ 1)  (44^+2(y^i^  +  a)^^) 
=  (44^  +  4a4^  -1)^+8  “ ®) (24  +  yct^+1  +  a) 


2  ( 3WJ. 

=  -= 


3  ; - j 


/ 


1-Kx' 


{Ja^+l+  24^X1+  44 Vo^+1  +  “  24^(3^?+1 -o+24^)(6'+^^^+24^ 


i^+l  -  a+24^)  (l+4^a^+l  4^+44^)  -  4( /x^+l-a)  (a+ Ja^+l+Zi^)  iJa^+'^+2^h '. 


^  +  e^(2+  4(o^+l)  -  2(3(a^+l)  -  +  2apZ.) 


+  4^(4ya^+l  +  sja^+l  -  4(4^1  +1)) 
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45. 


Co  f/  Calculation 
4 
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46. 


Let  h*  = 


★  * 

=  m  +  n 


(t+^c  ^u-t  “  uic^)  ’  tic  ■  uicO 

(-4  +  u“ic  4 

t^(u4.t)  t-iC  u+t  t+iC  u+t  J 


„*  1  /  -1  1  V 

”  “  ft  '  2  2  ~  ^^~~Z  ' 

®  t^(t+ic)-^  t^iC-t) 


Therefore  f 


^  (t^(t+iO^  t^(t-iC)^) 
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47. 


* 

=  p  -♦  q 


* 


)  -  2  (t(t+iO 


) 


>v 

f  = 


-1 


4(tt) 

1 


((t+iC)(M?))  ((t-iO(t+iC))‘ 


[g-An-Bg-KUA^-Bn)]^  ^  |  g  l.(l-fAg+BTi)  j 

4(a^+l)  l_[(A-q)2  +  (B+g)^J  ^  [(A+ti)^  +  (B+g)^]  ^ 


_ 1_ _ f  (g-An-Bg)^  -  (UA4-Bn)^  ^  (g+An-Bg)^  -  (UAg-fBn)^ 

((A+q)^  +  (B+g)2)2 


4(a^+l)  \((A-n)^  +  (B+g)2)^ 


^  i  ^  (g-An-Bg)(UAg-Bq)  ^  (g+Aq-Bg)  (UAg+Bq) 
2(g^+l)  V((A-q)^  +  (B+g)V  ((A+q)^  +  (B+g)2)2 

. _ 1 _  ^g^-1  -i-g(q^-g^)  +2qg-2g(A-foB)  -t-^gCB-oA) 

4^a  tl)  ^  (A-q)^  +  (B+O^r 


a^-l*a(q^-g^)>2qg-2g(A4gB)  -  2q(B~gA) 


((Atq)2  t  (B4g)2)2 


_ 1  [  2aT2(aA-B)g-2(gB4A)q-2aqg4q^-g^ 

4(g^4l)  \  ((A-q)^  4  (B4g)^)^ 

2  .2 


_  2a42(aA-B)  g42(gB4A)q42gqg4ri  -g 
((A4q)^  4  (B4g)^)^ 
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^  /  t  (t-lQ  ,  t\t+iO 


4(tt)^  l(t+iO(t-iO  (t-iO(t+iO 


1  (a-i)(a-ATi-Be-l(l+Ae-BTl))  ^  (g-t-l)  (a-^An-Bg-HKl+Ae^-Bnn 

4(a2+l)^/2  (A-t,)2  +  (B+e)2  (A+ti)^  +  (B+O^ 

j-.  (  a.^,-.V?(Al+BO-A4+BT]  ^  o^-l-Ki(ATi-Be)-Ae-Bq  j 

1)^/2  V  (A-n)^  +  (B+4)^  (A+ti)2  +  (B+e)2  / 


4(a^+l) 


-1 

4(a^+l) 


I  ATi-»B£-2a-o(Ag-Bn)  ATi-Bg+2o-4a(Ag+BTi)  ] 

^  V  +  <8+5)^  (A+l)^  +  (B+5)^  7 


1*  -  ~ 


A 

2tt  \  (t+ic)(t-ic)  (t-ic)(t+ic)  y 


t(triC) 


t(t+iC) 


-1 


1+Ag-Bn 

1  “1  .,2  + 


l4-Ag4^BTl 


\(A-’l)  +  (B+g)^  (A+n)^  +  (B+g)' 


g-ATi-Bg 


g->-AT-Bg 


A  =  B(  Ja^+1  +  g ) 
B  =  A(  J 0.^+1  ~  a) 
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* 


+ 


L _  [  a^-l+BT](l-a^-a</a^+l)  - 

^  (A-n)^  +  (B+e)^ 

a^-l-BTi(l-a^-a\/a^+l) -  AjiX-a^+l-HX'ja^+l)  | 
(A+B)^  +  (B+e)^  / 


i 

2  ,^3/2 


-2a+Ail(l-a^-Ki>|a%l)  +  Bg(l 


4(a  +1) 


(A-n)^  +  (B+^)^ 


\ 

^■KtJa^+l)  +  Bg(l~a^-a\/a^+l) 
(A+ti)^  +  (B+e)^ 


-2a“Ail(l-a 

.2  - .2 


2  , 

g  -  I  f  _ g 

2  ^  ,.r^2 


Bn+Ae 


A^-BB 


4(g  +1) 


4(g%l)  V(A-Tl)^  +  (B+0^  (A+B)^  +  (B+4)^ 


jg  {  B£  -  An  •¥  /g^+1  _  B£  +  An  +  Jct^+1 

^  4(g^+l)  V(A-n)^  +  (B+0^  (A+n)^  +  (B+0^ 


Let  r  = 


g^-1 


T  ‘3 
4(g  +1) 


f.  =  -  f 


3g^+  1 
^  ^  4(a^+l), 


*  *  g^-1  £ 

s  =  p  +  - 5 —  f„ 

4(g^+l) 


^  ^  ^ 

Therefore  =  f  +  r  +  s 
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Special  Cases 


1)  Re(f  )  « 


_1 _ (a^-l-Ht^-l+2a^-2  ^  -  l  +  o^  +  l-  2o 

4(0^1)  V  16(a^X)  _ 


o) 


A(a^+1)  \  4(a^+l) 


1  ^  +  o(  'Jahl  ♦  o) 

16((i^+l)^  8(a^+l) 


Re(s*) 


a 


^icJTTi-a) 

a(  Jok^+l  +  0-) 

8(a^+l) 


) 


Re(r*) 


4(a  +1)  I  ^2(a^l)  2  p7l  <.  J^l 

) 


-  o) 


) 


3a^+l 

8(a^+l)  Va^+1 


+  1  + 


/o^+l 


o(  Jct^+l  +  o)  _  7o^+l _ 3o^+l  _  o(3oi^-t-l) 


4(c?+l) 


16  (0^+1^  8(o^+l)  8(a^+l)^^^ 


0(1-0^) ,  J-  _ 1  3(11^-!) 

8(a+l7^  ®  4(a^+l)  16(ei^^l)^ 


_  l-g^  (  2o  ^  3  2  \  1 

“  CaW’a^iy'^ 
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ImiJ  )  = 


4(a^+l) 


( 


2a 


■fa  “  v/a^+l  f  a  f  /o^+l  +  a  +  a  +  o  + 


16(0+1) 


Itn(s*) 


2a 


8(a^+l) 


8<a^+l) 


4(a^+l) 


+  o  -  7a^+l  +  Q  -  \/o^+l  -.  o  -0+0  -  o 

/  2  ^ 

4(  yoSi  -  o) 

^ 

“  -  o  ^ 

+  /o^+l  +  ®  ^ 


0^+1 


2 


4  ,^0^ 


Vo^^+l  -  o 

K^o^+I  -  o)y 


Iin(r*) 


(a 

2  2 
8(0+1)^ 

3 


+  ^ 


+1) 


3o' 


4- 


a 


8(a^+l)^  8(o^+l)  4(o‘+l) 


2 

o 

2.,v3/2 


2)  Rein  =  - 


1  [  o^-l+OT)^  +  2n(B-oA)  o^»l  +  on^  -  2ti (B-M) 

4(0^1)  \  ((A-n)^  +  ((A+n)^  + 


((A-n)^  f  B^)(,A+n)^  +  B^)  -  B^  +  (A^-n^)^  +  2B^(A^+n^) 


2  4  2 

o  +  1  +  TJ  -  2oTl 


TECHNICAL  RESEARCH  GROUP 


52; 


401)11^  ♦  2il(ll-2Ami|(£^ 
<» /—2 . « \  /  -  2  «  ._A  «_^2%2 


2(a‘-+l)(o  +l+ni-2on  ) 


-(o^-l  +on^) (a^*l  42(2^^ -t-o^n^-^nS  +  n^) 

rrr7r7T7"i4TTiT  ^ 


2(a%l)  (o%l+n^-2on  ) 


2  ,  „2 

_ o  -l-H»n  _ 

2(o^+l)(a^+l+n^-2oin^) 


-  2 


"5 — — T" — ^ - TT 

(a  +1)  (o%l  +  TV -2oH^) 

_ 2£i2£t2^ 

2  ,  _4  «  _2v  TTTw72T77I5^ 


2(a^+l)  (0^+1+  Tf-2aTi^)  (o^+l)  (o*+l+ti^-2oTiV 


2  2 
d  +34on 

9  5  X  o 


(a^+l4a^-2oTi^)^  2(o%l)  (o%l+tr-2oti‘) 


Re(8*) 


aBTj 
4(o  +1) 


( — ^ 
\  (A-ti)^  + 


(A+ti) 


_  — 

4(a%l)(o  +l+Tr-2aTi^) 


m. 


2(o^+l)(a^+l+T]^-2aTl^) 


Re(r*) 


30L+1 


4(oi^+l)  tT^Zot^^-wi^I 
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53. 


fi.  -  5  T  2~? 

(a  +1+T1  -2aTl^)^ 


TT  2 

a  +1+TI  -2aTi 


1  /  x+£\  2  x+/ 

2  i  ■  / 


Itn(f*) 


4(a  +1) 
.2 


2  2 
2ot~2  2ot+2  (A+Bo) 'H  Hh.n 

((A-ti)^  +  ((A+n)^  +  B^)^ 


(2a-Hi^)4An(/o^+l  +n^)-2(A+Ba)n(a^+l+2(2\/a^+l  +a)Ti^.|.  nS 


2 (a^+l)(Ti  -2an^  +o^  +  1)  ^ 


_ ^(A+Bo)t] _ 

(a^+1)  (T]^-2ail^4a^+l) 

2knC^^  +  iJa^+l  +  2a)r\^  +  2a'Ja^+l)  -  4(A-t-Ba)n(ya^4-l  •4-a)Tl^ 
(a^+1)  (n^-2aTi^  +  a^+1)^ 

(A+Ba)(  /a^+l  +  o)  =  A(1-ks  yJa^+1  -  a^)(  \/o^+l  +a) 

=  A(  Ja^+1  +  2o) 

-(A+Ba)Ti _  2Ati(ti^-(  Ja^+l+2a)y\^+2a:  Ja^+l) 

2  4  2  2'*’  2  4  2  2  2 

(a  +1)  (T]^-2aT]'^-K»  +1)  (o  +1)  (T]^-2aTi  w  +1)^ 

_ (A-Bo)Ti  ^  2An(-n^^2or  Va^-fl) 
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54. 


Iui(r*) 


(3o^+l)  Ati(2o-  >la\l  -  A 

v/a^+l(Ti^-2aTi^  +0^+1) 
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55. 


3)  ReCf") 


1— 

2(a%l)  \  (A%(B+0  )  ^ 


(  ^  -  2B5)^-3^+l+4^+4B^e^+2  yo^+ie^-4B4(  y^?+l+  e^) 

=  a^+l+4\2e^(2  yja\l~a)  -  4B|(  v/a^+l+e^) 


Re(f*) 


_  -o)  4  +  UihJJ^+iha-ioJ^-oh 


2  4  2  2  5~ 

2(o  +l)(e  +2ar-wi  +1)^ 


^  2(A-KiB)e(o^+l+g\2e^(2£^ -  o))  -t-  (a^-l-oe^UBg(^^^^^^ 

2(0  +I)(§\2ar-«»>1)^ 


r*  r* 


2  2 

f*  .  o€  +  1-0 

®  2(a^+l)(?\2oe^4a^+l)  (o‘+l)(r+2ae‘-wi‘-+l)^ 


r.2  , w4  «  .2.  2  ,.2 


2  2 
06+1-0 
"5 — - 5 — T 


*  a.) 


J'  J 


2(0+1)  (6  +2a6  No  +1)  (o  +1)  (r+2cir-»o  +1)‘ 
2  2 

Qg^  -  3  -  o^  ^  2 

2  4  2  2'*'  A  2  5  2 

2(0  +l)(g  +2or-K»  +1)  (g  +2og^-«l  +1)^ 


Re (8*) 


oAg 


2(0^0)  I  A^+(B+g)^ 


oAg(Vo  +l  +  g‘--2Bg) 
2(o^+l)(g  +2og^-»o^+l) 
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2(a^+l)(e^+Le^-K»2+l) 


* 

"e  =  - 


3a^+l 


4(a‘’+l)  )  r+2ae  40  +1 


-4a^-4 


(e^+2ae^-K»^+l)^  2(a^+l)(4^+2a4^+a^+l) 


-11  2S±i\  x+i 

"  2  i  “  i 


f* _ Ag _ _ 

°  (a^+l)(e^+2ae^-Kt^+l)‘ 


(l-Hx/o^+l-a^ClVaiWa^+l  -  0)40^4-1) 


-2(Ja^+l  -a) (a|^4-(a Ja^+1  4-  1-a^) ^^+(1-0^)  Ja^+1)^ 


AgCl4o^-Q^^4-l)  2(l-a^-ictJahl)  Ja^+liK(aKl)(l 


(ot  +l)(4^+2ci^  ■+Cfc  +1) 


-  0^4-0  v4^4  ])  +2(0^-  lXa^4-l-o 


'J 0,^+1  -  o) 


fo+l  (  ^^20  5^'4a^4-l) 


Ag(2  ■  2(o^4-l^ 

(0^4-1)  (1^4-20 1^40^4-1)^ 


+io^4-l) 


°  2(a^4-l)(e^4-2oe^4a^4-l) 


,*  -<3oSl}  * 

*  ^  2A|<;i^2a->/Si; 

^  /  O  /.  O  /> 


A£(ag^-oio^4-l4-2(a^4-l^ 


^  2(071) (|\2a|2^^1)  4(o7l) 


3a  4-1  ,* 
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57. 


4)  -4(o^+l)f*  -  -27T^ 


((o^-l-K>(n^-e^)-2|(A4aB))(D‘(tl)+D{’(-n)) 


D{(n)D‘(-n) 


+  2n(4+B-aA)(D^(Ti)-D^(-n))) 


where 

dJ(ti)  =  (aW+(B+0^2Ati)2 

1  (A^+  t?+(B+0  ^  +  AA^n^+^An  (A^+  ^ +(B+0 

2(a^+l)f*  -  -:5 - h -  ((2oi^-l-2e(A-KiB))((N£^+o+2e^+2BD^ 

D^n)Di(-n) 

+  2C^?^  +  a)(l^-K»))  +  8A(e^-K»)(e+B-oA)(^^4o+2e^+2B|)) 

D  (il)D/-n)  -  (>^!^+<»+2e^+2Be+2An)(^^  +  a+2e^+2B 

-  (£^+a+2e^+2BO^  -  2(^i^+ci)U^+o) 

-  4^\8Bt^+A^a^+l  +  4(^o^+l  +oi)BS+l 

2  2o^-1’2£(A-k»B1  ^  — 4  - ((  2o^-1-2|(A-»oB)) (^0^+140) (i^-K>) 

-2(o+l)f  . 

4  Cjo^+l  +o+2e^+2Be)  (2A|+l-o(ioi^+l  +  »))) 

j 

^  2o^-l-2fe(A-K»B)  ^  (^2+1+  o)  (o(o-\t^+l)-2oB0 

D^(Tl)D^(-n)  Dj(T|)Dj(-n) 

+  2(e^-Ki)(|^+Bl5)(2Ae+l-a(Jo^+l  +«))) 
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Dl(n)D3^(-Tl)  Di(n)D^(-T)) 

+  2(2A$^+(2-o(^!^  +  o))e^+(l-o(^i^+o))BO 


2o^-l-2fe(A-K»B)  ^ 

Dj^(n)Dj^(-Tl) 

2a^-l-2fe(A4oB)  , 
D^(n)D^(-Tj) 

+ 

D^(Ti)D^(-n) 


(45^+2BC^+(£^.4a)(4Be 

D^(n)D‘(-n) 

T ■  I  (4e\2Be^+3oie^+2Bae-a^) 
D‘(tOD‘(-tO 

(4B€|^-2cie^-2o^e-^^B) 


^-2ae2,^«)) 


2o^~l-2g(A-KiB)-f4-f4(yo^->-l-Mi)Bg 

Dj^(Ti)Dj^(-n) 

+  ~5 - -  (-6Be\(3oi-4,f?^)e^-  (4^?+1  +  2o)B§-(o^+1)) 

P£(n)D^(-Ti) 

^4jEi±21( (2a-4^^) e^-4B^^  e^-2(o^+l) e- (i+4a^)B) 
D‘(ti)D‘(-t,) 

2o^-l-2g(A-K»B)-t-4-»4(Xi^-fl-Ki)Bg 

Di<n)D^(-n) 

-  —5 — -  (4^\6B |\(4i*^+l-3oi)  {^+(4^^+l+2o)B |+(a^+l)^ 

Di(n)Di<-Tl)  ' 

-  (2oe^+4B^^  t^+2(o^+l)4^+(^^+l)BO 

D‘(Ti)D“(-n) 
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59. 


-2(a^+l)f* 


Let 

•-2(a^+l)f*  = 


2o^-l“2g(A-K>B)-t-4(s/a^+l  +a)Bi~2aiJQ^+l-¥  a) 

Di(n)Dj^(-n) 

-  -5 - ^ -  (-2B4^-3oe^-2aB§4a^) 

D^(Tl)D‘(-Ti) 

-  CUB(Ja^+l  -  a)  ihlJa^+uJoKl-a) 

D^Cn)  D£(-ti) 

+  (l+2a(a-^i^))Be-|  ) 

-l-2o^/oVl■^^2^iL^+l  Be 
D^(Tl)D^(-Tl) 

+  -X - ^ -  (2B5^+(2^!^-3oi)e^+  (L^-3o)Be-  ft(£^-a)) 

D‘(ti)DJ(-ti)  2  V 

Di(t1)Dj^(-ti)  -  D2(e) 


2^^+l 


ir^- 


(4Br+2(2( 


n^iOD-i-o 

+  2((Ja^+l~a)  -2a)B4-  (^+l-a)o)D,  (-^ 


JaKl-a)i" 


^Y°»  +^.LBS  q;  ^ .  -y— ^ - (2Br-oe-2oB) 

D2<0  DJ(0D2(-0 

.  ■ -  (vf^-o) (4e2+2B§-o) 
D£(0D2(-?) 

D2(-e)  -  D2(0-  16Be^-8(^^+o)Be 
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60. 


■2(a2+l)f*  - 


.  4£(2Be  -ag-2aB) 

(O  D2(0D2(“^) 


^  (2B|2-a5.2aB) 

D2(0D2(-0 


.  2(>/a^+l>a)(4£^-h2Bg-o)  ^ 

*  f\  /ASrT  /  ifc\  +  O 


Do(ODo(-0 


1)2(002(0 


Nj^  5  324^(2§^+('^i^4a))((>£  ^+l-a)  ^^-oB^-a  (>/o^+l-a)  ) 

-  324^ (2(^^+l-a)  4^-2aBe^+(l-2a (\/o^+l-a) )  ( /a^+l-K»)B|-a) 

«  16§^(\/a^4l-a)  (44^“4(^a^+l-Mi)oB4^+2(yo^+l-o)4^ 
-2a(^^fl4a)^B|-2a(^a^+l4a)) 

N2  a  16Be(^i^-a)(2e^+£^4<»)(4e^+2B|-ci) 


I6B5  (/o*"  f  1-0)  (84%4B§"'+(-2af4(\/o^+l+o)  )  ^^+2 


.3 


(  io^4-l4o)  B|  -a  (  /o^+l-K»)  ) 


Nj^+  =  16(7a^-l-a)  (e^(4e\8B§^-4a(>^?U4a)B|^+4>jKie^-4a|^ 

+  4(\/o^4l4a)B^-2o(l+(^o^+l-wi)^)B5+l-2o  (70^+140)) -B4o(/a^+l-Wi)) 

-  16(^i^-o)(e^D2(0-4oe^-2oBO 

4  (7a^+l+a)  (§^(-4aB|^-2oi(^^4l4a)B§-2a)-aBO) 

=  16((^Kl-a)  (5^02(0 -^?^-2aB5^)-ae(4B§^4-2B(^i^4a)  4^4244-8)) 
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-ae(BD2(0  -4(yoi^+l-a)  4^+2(a- 

-  16(((,PTI-a)4^-oBOD2<0) 

_2(a2+i)f*  , 

©2(0 

+  -4£(2Bfe^-qe-2oB)-2(ia^-t-l-o)(4^^-t-2B£-o)-»-16((jti^-fl-o)t^-oB6) 

D2(e)D2(-e) 

-2(a^+l) £*02(0^2  (-0  -  (44^-8B4^+4jKIe^-4(^^4a)Be+l)(2^i^B|-2o^o?a-] 

-  8B4\4o?^+8(^^+l-o)e^-4(>/o^+l-Ki)Be+2o(i:i^+l-a) 

Let  Q  -/ 4§^  •♦: 

-2(a^+l)fV  -  Q(2j!»^+lB4-2o/a^+l-l)+2^^+l(-4(ia^+l-a)e^ 

+  4(>/a^+l-o)Be^-2e^+2BO+16o/o^+lBe^ 

+  4^[^(1  -2a<^i^+l-o)  )  ^^+8o\^i^+l(  Ja^^-Ki)B5 
-  2 Ja^+l (ofc+i/o^+l) 

-  Q(2io^+lBe-l-2(o^+l))+8\/oi^+l(\/o^+l-K»)BC^ 

+  4^^+l(l+2a  (\^^+l-K»)  )B|-4(a^+l) 
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-2(a^+l)fV  -  Q(2\/o^+lBe-l-2(o^+l))+8/o^+lA5^ 
+  4vj£i^+l(ia^+l-Ki)A|-4(o^+l) 


Maliii  * 


a 


a?(Dj^(ti)  +  D^(-n))+  bti(d^(t))  -  D^(-n)) 

0,(0 


liafiii  *  _  .  A^(Xi^+l-wn-2g^+2Bg)+g' 

o  "  ^2(0 


40 


2Ae^2eV^ 

D2( 


L|yQ^4.1-fOl)  fe40 

g*Q2  .  <2Ae^+2^^+A(^^^+l-^o)  |4o)  (Q-SBe^+Acio^+l-o)  |^-4A^+2) 

-  (2A5^+2e^+A(^i^4o)  e^o)Q-8e®+8  (  (^io^-o) A-2B)  4^ 

+  (8(^i^-4»)-8((>^?+iL4o))5^+(4A-8A+4A-8oB)5^ 

+  (444o(^^+l“o)-2(/o^+l4a)^)5^ 

+  (2(Ja^+l'Ki)  “^)A§+2o 

-  (2A5^+2e^+A(^?+l4o)  e-Wi)Q-8e®-8B^-16o4^-8oBe^ 

+  (2-8a^)e^+2(io^+l-  a)A^4-2a 
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g*Q2  _  (2Ae^+(A(^?^4o)-2B>5-Ki)Q-8oe^-8oi^5^+(2(^^-a)A-2B)e+2o 
-  (2A|\(3o-/a^+l)Ae-a)Q 

\ 

2(o^+l)(f*+ s*)Q^  -  (2Aae^+((3a-^^+l)oA-2^?^)e-a^+l+2(o^+l))Q 

r  8  JL^+U4^-4^^+1(Jo^+1-ki) Ae+4(a^+l) 

-  (2Aoie^+(a -2) Ae4a^+3)  Q-s/i^+UI^ 

*i-  4 ya^+l(j[»^+l-wi)A444(oi^+l) 

2(a^+l)r*  -  (3a^+l)  ^  (1  -  -ML  -  (i  +  -~z“  > 

yja^+1  Ja^ 

2(a^+l)f.Q^  -  (2(a--^^^^  Ae^+  (a(/i^^a)-3(a^+l)-  °*^j^^)Ae 

+  4a^+4)Q  -  uJahl(2Ai\iJaKl-M)kO-¥UiaKl) 

-  (2(a£^-3ci^-l)  -  (2oi^-»^^i^(2a^3)) 

+  4(a^+l))Q-4^^+l(2Ae^+(yo^+l-»o)Ae)+4(o^+l) 
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appendix  II  -  J^(a) 


Let 


Let 


•  2, 

u  du 


-oo  P 


2  2 

where  P(u,a)  ■  1  +  (u  +o) 


|i  =  P‘  =  4u(u^  +  a) 
Integrate  by  parts 


1  -a  !  — ^  p'  da 

"  .«  p"+l 


4„  /  “i<»L±gl  du 

-  p"+l 


2/  2  .  ■  2  ,  2 
u  (u  +  a)  «  P  -  au  -1-0 


Therefore 


or 


h+1 

2^ 


h+1 
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Finally 


n 


-T/ 

-00  Jf 


4  2  2 

-  P  -  2au  -a  -  1 


Therefore 


3J 


«*  ..2^  2. 


p 


since 


«  2,  2  .  I 

/  “  <“>)  du  -7^1 

L  p”+i  *” 


and 


therefore 


(l«^)I„^l  -  (1  -  i)l„  - 


o  al 

(1  -  7^)J  -  J  .  +  -7-S 

'  4n'  n  n+1  4n 


For  recursion  purposes 

al 


Vl  -  •  "SF  +  -  ^J-'n 


‘•n+1 


1+a 


S  -  b\  -  “J„+i) 
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n  »  1  Determination 


P  ■  (u-z)  (u-z)  (u+z)  (u+a) 


Where  z  and  -z  are  in  the  upper  half  plane,  and  z 


7  ^  .  0 

J  p 


Therefore  Jj^-zzl^^-O 
i.e.  «  v/a^+1 


du 


.  2  2.  .  2  -2, 

(u  -z  )(u  -z^ 

1<  ■  2Ti  (  2  «.2  "  ^  _2  2  ^ 

2z(z  -z  ;  2z(z^-z^) 


Since  z^  -  z^  -  2i 


Ii(a) 


2  zz 
irB 


I''*'  . 

/a^+l 


irA 


■  -a+i 
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APPENDIX  III  -  Steady  State  Integrals 

Let  L  ,,<o)  -  /  2 -  r  '  dx 

where  P  >  F(x,a)  Is  defined  in  Appendix  II. 

We  wish  to  evaluate  and  l,2» 

Using  Appendix  II  and  suppressing  the  argtiment  o  In  1^^,  J^. 


*^2  ^^2 


J2  + 


a 


31, 


1-KS 


I  (—  -  aJ2) 


3al, 


140 

-al 


4(140^) 


4(l4a")  4(l4a‘') 

2a  +  \laKl  T 

4(l4a2) 


ttB 


4(140"-) 

I, 


(1  + 


2a 


3oI, 


4(140^) 


4-1 


ttB 

4^140^ 
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appendix  IV  -  Condition  (4)  Integrals 


Let  P(|,a)  be  defined  as  in  appendix  II. 


Let  M 


c,m 


;  —cos  (— )  d? 


.oo  P 


,  (-1)”  ,^,2n  r“  de 

^  TW  J  ;si+i 

o'  '  “OO  P 

I  (-i)"(2n)^"  , 

(2n) !  2n+in 

o 


Let 


20 

N  "  /  «  cos  (*=-) 


“  (2n)i  '^2n+iD 


Let 


1  2o 

M  =  /  —  sin  (^) 

s.m  •*  „ni  'P'  ^ 

-  go  r 


Let 


I  I 

^  (2n+l)l  ^2n+l-Hn 

o 


-s.™  ■  /  ^  ->« 


”  i^iAzsa^ 

^  r9nxlM  *^5 


(2n+l)!  ‘'2n+Uin 
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r»  cos  — 

l,c  =  P 


di 


-  ^  <2«)"”<-'2n.2«hn+2) 

O'  ' 

I  /  i!is^  CO.  f  de 


1  ;  (-l>”(20)^”  ^2ntl 

%  ^  <inr!  2o+T 

O 


Sltallarly 


1  2  ^2n+2 

o 


'2,8  “  4  t  (2n4l}!  2n+2 


■t..  ■  “•  T  '« 


1  r  d  / 1  2iiv  gy  j* 

l“c.2-^“  ““  d«^  * 


2ftv  dP 
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Since  P  ■  (§  )  +  2a(^  )  +  a  +  1 


^  +  2a 


din 


Therefore 

* 


■l.c 


d  /  1  2n.  2  2n  ^  2n  2n 

^  CO.  — )  -  .  CO.  ^  .in  ^ 


1  M,,, .  r  ^  (S  .in  f  -  c  f)  de 


-  T  “c.2  *  *«  t“<''..4+'>“.,4>-<''c,3  +  »“c,3>l 


n/«^v2n 


- 1  * 


^2n+2 


2n<»2 

2n-»'2 


+  W2n+5«hn+5> 

o 


»n/««\2n 


*I,c  -  *  ‘*"+«*2»42> 
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*  i 
^1,8  “  2 


1 

“  2 


o 


QO 

=  z 

o 


r  4  2fi  j* 

si" 


M  -  +  o  r  4  (4  sin  fi)(2e^  +  2a)<l{ 

8j4  p 

M  2  •  ^  (  I  cos  T 

8)4  _«o 


“5,2  -'*“‘"^''c,4+“*'c.A>  +“.,3+““.,3’ 


00 

z 


.  rl  ^  r,  ,  T  >.ls±l 

•  ^®L2  ^  {airr  ^'^2n+4^“’-2n+4^  2n+l 


-  ^<"  -  |)«2«4-  “h«A» 


/  i\n/rtrt\2n+l 

2  il,c 

1  ^1,8 
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Let  E  ^  2ci  */oi  4*1 


.  r  (M) 

J7Z^  -- 


2A_  “  (-i)”(2n)^”  rliiif+aiiiisi  di 
7-  \  (2”)>  i»  V"+5  ' 


2A  *  (-n"(2n)^"  .^2n+2  ^^2ih-2  v 

"  ~7^  ^  2n!  ^4(2n+2  ■•■  4(2n+2)^ 

^  °  n- 

A  -  (-l)"(2ft)^”  /^2n+2  +  (2ojrOjfl)I 

‘  7^  f  (2n)!  Z55+2 

2^0  +1 


2n-f2 


Slmilaxly 
1 


3»®  r2T7  o 


2/a>l 


^  12^ 


3J,_,+ (2a-yKl)I 


2n4.3"  V—  ) 


2n+3 


Let  F  ■  a  ja^+1-  2(o^+l) 


.  _2A^  j.”  CO.  (f  )  dt 

-*  * 

.2, .2 


^  CO.  (f)dt- Ij 


2(a%l)  -•» 


4(0+1) 


-  +  U2  +  U3 
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T-  t  (2ni!  {„  p2''+*  ' 


2A  "  .^2n+3  ^  ^^2n-t-3  x 

^  ^  i  (2n)!  Wn+3)  4(2n+3)^ 

/a%l  ° 

_A  I  oe 

2(a%l)  o 


-  OQ 


A  ”  (-l)”(2n)^”  .^2n+2  ^i2s+2_^ 

I  '  'inf'  <4(2n+2)  ‘  4(2ii+i>^ 

8(„i+i;^  ^  2n+2  ^ 


u. 


8(a%l)'’'^  o 
-  (f/o^+I  +  (3ci^+l)E)l2^2^ 


2ii+2 


'•211+2' 

+1  +  (3o^+l)E  -a(o^+l)  -  2(oi^+l)^^+  (3o^+l)(2a- 


.  o(7a*+3)  -/b^+1  (5o*+  3) 


.i.,>J/4  ^  (2n>! 


(-I>'*(2a0^”  ,3tJahl  -  (3a^-H)a., 

^  2n+2  ^ 


8(a  +1) 

(5g^+^^^+l  -  o(7o^+3)  , 

+  ,  O  ^ 


’2n+2 


2n+2 


2n+2 


'211+3 


2n+3' 
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Similarly 


n+1)!  ' ln+3 


(50^4-3)  Jb^+l«o  (70^+31 

2n+3 


'2n+3 


+ 


4(o^+l) 

2n+4 


^^2n+4 


+ 


(Ta-Ja^n^^)) 
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APPENDIX  V  -  Force  and  Mmant  Terms 
Parc  1. 

"Sc  -  /  t  ■*’' 

«  00 

•  00 

®k  ■  /  + 

M  09 

2  2 

Let  P  ■  P(Ti,-a)  -  1  +  (ti  -a)  (See  Appendix  II) 

Ri  -  8i  J  d, 

»  81  (J2(-a)-al2(-o)) 

Henceforth  we  will  suppress  -o  in  and  I^ 


-al2  -  -=^  (r^  +  0J2) 

^  1-KX^  ^  ^ 


31, 


(1- 


1-wi 


2")J^ 


1-KX 


1  *^1 


R, 


21  ( 


2a 


1-w 


2)Ii 
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*1  ■  ^ 

*  -fl»  TJ  P 

-  -'“<"2  ■  “  L  ’ 

/  ■  ■  n?  ^  ^ 


dn 


nP 

•o  2 

-8  f  3^  dn 

ioo  p-^ 


.  -803-013) 

R*  .  -4je(i2  +  8003  -  0I3)) 


0X3  -  (0J3  +  ^  I2) 

1-K» 

-  I  "  .  2^8  ”*■  8  ^ 

140  l-to 


R*  -  -4i(I«  +  (5J2  -  60I2)) 

^  2  140^  ^ 

1-50^  (oJ  4-—) 
XT  +  4  ^ 


(1-  -^)I. 


140 


(l-»o‘) 


-SL-  (5  +  1=^)J2  -  — 

140^  140^  /  ^ 


R*  =  -3i  ( - 2|-372  +  -2)11 

^  (1-K»  ) 
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T  2,  2  . 

.8i  /  2Li3^  dr, 

-rf  p‘‘ 


-2il, 


-^R 

2 


1-K* 


(2a- 


16Ai 


2,^2 


/ 


/  n  -«^0  -gl  dr, 


a+l 


-  *9 


n^exe  E  «  2a  -  Ja  +1 


2ML.  (1  J  .  E  I  ) 

-  4  ‘*2  8  V 


Jahl 


- ^  (0J2  +  I  ip 


8(140  ) 


aE 


8(140^) 


8(140*’) 


^  3  ((3(140^)  -  E)J^  4-  (3o(l40^)-(o^4.3)E)ip 


2(l4o2)' 

(oJ^4-  (0^4-  3)ipE  -  (o/o^4-l  4-0^4  3)(2a-^^4-l)Ij^ 

-  (2o (0^4-3) -a  (0^4-1)+ JKi  (0^-3)) I 

■  (o(a^4-5)  4-  Ja^+1  (a^-3))Ij^ 

3(140^)  (Jj^4oip  -  3(140^)  iJa^+1 4- a)I 
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»3  - 


(o  (20^-2)  +  1^1  (2a2+6)) 

i{x-wry'^ 


.  (o(o^-l)  +  ^0^+1  (a^+3)) 

- -  (o^(o^-l)  +  (a^+l)(a^+3)  +  o\4i^+l(2oi^+2)) 


2(140^)^ 


(2a^  +  3o^  +  3  +  2a(o^+l)^^^) 


2(Ho‘) 


R 


4 


l6Ai  f  TL^ni-g^^iE}  ^ 


where  F 


aJP+1- 


r"  n"(n%Kai^F)  .  _3c^ 

(a  +1)  P"*  4(0  +1)  ■* 


2(1-ki^) 


16Ai  ^3  ®^3v  4Ai  .^*^2  ^2x  3o^-i-l 


EI3 

IT 


b  +1 

Tl  +  +  I3 

— S-j^  (0J3  +  I  Ij) 


12(140-) 

oE 


12(140'^)  ^  12(1^0^)  **  ® 
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6(1+0 


-1-172  (5(3(1-k>^)  -  aE)J2+  (3a(lw  )  -  E(a  +7))l2) 


T, +T,  -  - ((9a^Ji^+ 5(3(140^)  -  oE))J. 


+  (3(/l4a‘  F+a(l4a^))-  E(o^+7))l2) 

9aia^+l+  5(3(l4a^)  -  oE)  -  l^<jck^+l+  5(a^+3) 

I;;(a^+7)=--3(2«(l-i<i^)  -  2(1h<i^') ^^^)  +  (ia^+l -  2a)(a^+7) 


2a(2a^-4)  +  /o^+l  (l-5a^) 


(2a(2a^-4)  +  /a^+l(l-5a^))l2  -  (oj^  +  _^) 

0+1 


9a^+  12q^+  15  +o^c?+r(9o^+  15)  (°h  ^  j 


0*^+1 


T1  +  T2 


8(140 


((3o^  +  4o^  +  5  +  o/o^+l  (3o  +5))J2^ 


+  (a(3a^  +  80^  -3)  +  /l^  (3ei^  +  l))Ii) 

Aili  A  9  /  9  A  9 

■■  i^w  (a(6o^  +  16o^  +  2)  +  70  +1  (6a  +40  +6)) 

8(1400^'^ 
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-(3a  +1)  All,  2  rr~  2  ,v 

T,  -  - *  (a(2o  -  2)  +  /a  +1  (2a  +6)) 

3  8(l-»a^)^'^ 


-All 


8(140  ) 
All 


^  gy2  (a(6a^  -  -  2)  +  /a^+1  (6a  +  20a  +6)) 


^  8(140^) 


(a(20a^  +  4)  - 


vi(Ja^+l±^  (16o^(a  -  Ja^+1)  +  4a(o^4-l)) 
16(140'^)^ 

- 2l_  (a(a^4-l)(v/o^4-l  4-  a)  -  4o^) 

4(140'^)'* 


R  -  I  c 


ko\ 


aioBI 


C  ,  rtR 


10*"!  - 


i  (2a^?+l  4-  1  -  3a^)(Ja^+l  -  a) 


0  4-1 
31oBI 


0  4-1 
aBlI 


i  (oi(5o^4-1)  4-  (1  -  5o^)) 


^20^2  “ 

aBll 


-  (2a  -  ^0^4-1  (o/o^4-l  4-  1  -  a^)) 


2 (o(l  -  3a^)  4-  \/o^4-l  (3o^-  1)) 


a‘-4- 1 
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aBn, 


R  =  -5 — -  (a (12a  +4)  +  /a  +1  (2-  12a  )) 
“  a^  + 1 


(«(6a^+2)  +  U^+1  (1-  6a^)) 

(a  +1) 


- ^  ^^3/2  +a)^((l-  6a^)  ( vjEt^+1  -  a)  +  3a) 

a  (a  +1) 


u MjgTT  //TT.  V,,  .2  rr 

(ya  +l  +  a)(l-  3a  +  3ava  h 


a(a^+l)^^^ 


-U^MiTT 

a(a2.1)3/2 


(4a  +  Ja^+l) 


h  =  J 


-  <30  2 

2  r  r\-a 


dn 


-.to  p- 


■al3  = 


=  16i  (J3  -  al3) 
-a  J 


140 
2 


2  (g  I2  +  0,3^) 


al. 


(1-  -2— )J  -  ^  j  + 
140  I40 


2/ 


“  ,  2  ^^*^2  ■  ^^^2) 


l40 

31 


ll^Jr 


Rn  - 


2(1-ki^)  ^  l4o^ 


3r  ,  1  2a  X,  r  ,  rrr 

(  “  - “  2^ 


-/i,  (  ' 


2(l4o‘) 
a  3 


/  1  .  ct  3  6a 

^  2(140^)  (140^)^/^  (140^)^ 
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S*  =  -8i^(l3  +  12a (J^  -  al^)) 


T  -g  /II  T  T  \ 

■“^4  ",  2  ^12  ’■3  ^  “■^4^ 


1-K» 


l4g 
2 

P"4 


g  X,  1  /“’■a  .  3  ,  . 

(.1  "  9)*^/.  “  2  ^12  ^4  *^3^ 


1+a 


S*  =  -8^^(l3  +  (9J3  -  IO0I3)) 

l4a 


•9a 


8u 


2(1-ki^) 


2 


l40 


I3  =  ^^2  <8  I2  “J3) 


_ L-  (740  ,  +  ^  J  ) 

1^2  ^  8  ^2^8  ^^2^ 

2 

(740^)12  =  (|  +  aJ2) 

140 
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